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Abstract 

We describe a parametrization of the real spectral locus of the two- 
parametric family of PT-symmetric quartic oscillators. For this family, 
we find a parameter region where all eigenvalues are real, extending 
the results of Dorey, Dunning, Tateo and Shin. 
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1 A family of quartic oscillators 

We consider the eigenvalue problem in the complex plane 

w" + (C^ + 2hC + 2iJC + X)w = 0, ^ 0, t +00. (1) 

Here J and b are parameters. This two-parametric family is interesting for 
several reasons. 

When 2 J is an integer, 2J < 1, and b > 0, the problem has the same 
spectrum as a spherically symmetric quartic oscillator in R'^. In this case 
2J = 2 — 21 — d, where I is the angular momentum quantum number [5, 4]. 

When J is a positive integer, problem (1) is quasi-exactly solvable (QES) 
[2]. This means that there are J eigenfunctions of the form 

t/7(C)=p(C)exp(-<V3-z6C), 
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where p is a polynomial of degree J — 1 in C whose coefficients are algebraic 
functions in b. 

When J and b are real, the problem is PT-symmetric. The eigenvalues 
of a PT-symmetric problem can be either real or come in complex conjugate 
pairs. Both possibilities can be present for J > 1. A very interesting feature 
is level crossing in the real domain: for some real b and J the graphs of the 
eigenvalues \k{b) can be real and cross each other. This phenomenon was 
discovered by Bender and Boettcher [2] numerically, then it was studied in 
[12] , where the presence of infinitely many such real level crossing points was 
proved for positive odd J. 

When J — > +00, the QES part of the spectral locus approximates the 
whole spectral locus of the PT-symmetric cubic family 

-w" + (i/ + iaz)w = Aw, w{±oo) = 0, (2) 

which was subject of intensive research, see, for example, [3, 7, 8, 11, 16, 11, 
20, 23]. 

By the change of the independent variable z — iC, problem (1) is equivalent 

to 

L,,j{y) = -y"+{z^-2bz^ + 2Jz)y = Xy, |/(te±-^/3) ^ 0, t ^ -hoo. (3) 

Shin's theorem [20] applies to these eigenvalue problems when J < 0, and 
implies that for J < all eigenvalues are real. The proof of Shin's theorem 
is based on the remarkable ODE-IM correspondence of Dorey, Dunning and 
Tateo [8]. Here we extend this result of Shin. 

Theorem 1. All eigenvalues of (1) or (3) are real for J <1. 

The condition J < 1 is exact, because it is known that for every J > 1 
there are non-real eigenvalues [12] . Our proof of Theorem 1 is based on purely 
topological arguments. Using the formulation (3), we establish a certain 
property of eigenf unctions for J = 0, and then show that this property 
persists for J < 1 and prevents level crossing. 

The real spectral locus ^(R) C is defined as the set of all real triples 
(6, J, A) for which there exists y ^ satisfying (3). This is a non-singular 
analytic surface in R^. The main result of this paper is a parametrization of a 
part of ^(R) corresponding to integer J in terms of Nevanlinna parameters. 
In [11] we obtained similar parametrization of the real spectral locus of (2) 
and another two-parametric family of quartics. 
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The family considered in this paper is much more comphcated because 
of the presence of QES part and of the real level crossings. In [13] we 
parametrized the real quasi-exactly solvable locus Z^^^{R) of (3) which con- 
sists of all triples (6, J, A) G R'^ for which there exists a function y{z) = 
p{z) exp(2;^/3 — bz) satisfying (3) with a polynomial p. 

The plan of the paper is the following. In the next section we introduce 
Nevanlinna parameters and state our principal result, Theorem 2, about the 
correspondence between parameters {b, J, A) G ^(R) and Nevanlinna param- 
eters. In section 3 we prove the "easy", algebraic part of this correspondence. 
In section 4 we study the case J — and in section 6 we prove Theorem 1. 
Then in section 7 we describe the parametrization of the part of the real 
spectral locus ^(R) where J is an integer in terms of Nevanlinna parameters 
and complete the proof of Theorem 2. 

In a forthcoming paper we will parametrize the whole two-dimensional 
real spectral locus of the family (3) in terms of Nevanlinna parameters. 

2 Nevanlinna parameters 

The references for this section are [18, 21, 9, 10, 11]. 

Suppose that b, J and A are real. Let y be an eigenfunction of (3). Then 
y*{z) = y{z) satisfies the same differential equation and the same boundary 
conditions (3), so y* = cy. We can choose y so that y{xo) G R\{0} for some 
real xo- Substituting this xq to y* — cy we conclude that c — 1. So our 
eigenfunction y is real. It is defined up to multiplication by a real constant. 

Let yi be a solution of the differential equation in (3), which is real and 
linearly independent of y. Consider the meromorphic function f = y/yi- It 
is a Nevanlinna function, which means that / has no critical points in C and 
the only singularities of are finitely many logarithmic branch points. 

In the case that yi is normalized by yi{x) — >■ 0, a: — )■ -|-oo, a; G R, we call 
/ a normalized Nevanlinna function of (3). The normalized Nevanlinna func- 
tion is defined up to a real multiple. Existence of yi with such normalization 
is guaranteed by a theorem of Sibuya [21]. 

Different Nevanlinna functions associated with the same point {b, J, A) G 
Z(R) are related by a real fractional-linear transformation of the form / i— >■ 
af/{f — P) where a ^ and /3 are real. 

Nevanlinna functions / of (3) have no critical points, because /' = {y'yi — 
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yy'i)/yl and y'yi — yy[ — const. They have 6 asymptotic values in the sectors 

Sj = {te'^ : t > 0, 1^ - 7rj/3| < 7r/6}, j = 0, . . . , 5. 

In what follows j is understood as a residue modulo 6. These sectors are 
in one-to-one correspondence with logarithmic branch points of /^^. The 
asymptotic values of the normalized Nevanlinna function are in 5*1 and 
S-i, because of the boundary condition, and oo in 5*0 because of the normal- 
ization of yi. We denote by c and a the asymptotic values of the normalized 
Nevanlinna function in 5*2 and S3, respectively. As / is real, the asymptotic 
value in S'_2 is c, and a is real. 

These asymptotic values a, c are called the Nevanlinna parameters. They 
are related to the Stokes multipliers by simple formulas [21, 17]. Normalized 
function / and the Nevanlinna parameters are defined modulo multiplication 
by a real non-zero number, so we can further normalize them. We will use 
different normalizations, depending on the situation. 

Nevanlinna parameters (c, a) e C x (RU{oo}) are subject to the following 
conditions: 

c^O, c^a, (4) 

and the set {0, 00, a, c,c} contains at least 3 distinct points, in other words, 
the combination c = 00, a = is prohibited [18, 21]. 

Nevanlinna parameters (modulo multiplication by real constants) serve as 
local coordinates on the real spectral locus Z(R.), [1, 17]. Notice that the set 
of pairs satisfying (4) modulo proportionality is a non-Hausdorff manifold. 
But this will cause no difficulties as we always work in local charts. 

Relation between (6, J, A) and (a, c) is very complicated: Nevanlinna's 
construction of the map (a, c) — )■ (6, J, A) involves the uniformization theo- 
rem. So it is interesting and challenging to establish any explicit correspon- 
dences between the sets in the space of Nevanlinna parameters and the sets 
in the (6, J, A) space. 

Our main result in this direction is 

Theorem 2. J is an integer if and only if 

a — or c — c. 

In the next section we prove the "only if" part. 
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3 QES locus and Darboux transform 

In this section we prove the "easy part" of Theorem 2: if J is an integer, 
then either a — Q or c — c. 

Suppose that for some (6, J, A) e we have a = 0. This means 

that the cigcnfunction y tends to zero in 5*3 (and also in 5*1, S-i), so y is an 
elementary function of the form pe^ with polynomials p and q. Substitution 
to (3) gives 

y{z) = p{z) exp{z^/3 - bz), 

where p is a polynomial. It is known that such eigenf unctions exist if and 
only if J is a positive integer [2, 12]. Such points (6, J, A) form the QES 
spectral locus Z^^^(R) which consists of smooth algebraic curves 

Qj{b,X) = 0, J = n+1, n>0, 

where Qj are real polynomials of degree J in A. QES spectral locus was 
studied in [12, 13]; in the second paper it was parametrized in terms of 
Nevanlinna parameters. The converse is evident: if (6, J, A) e Z^^^CR), 
then a = 0. 

In [12] we obtained the following results: 

(i) If J is a positive integer, and b is real, then all non-QES eigenvalues are 
real. 

Let Zj(R) = {(6, A) e R2 : (6, J, A) G Z(R)} and let Z^^'^(R) be simi- 
larly defined. Let Z*j be the closure of Zj(R)\Z^^^(R) in R^. 

(ii) When J is even, Z}(R) nZ^^'^(R) = 0. When J is odd, then c = c holds 
at all points {b, A) of this intersection. 

(Hi) If J is a positive integer, then Zj = Z^j. 

To deal with the condition c = c, we use the Darboux transform [6, 14, 12], 
which we recall. Let ipg, . . . ,ipnhe some eigenfunctions of a differential oper- 
ator L — —d'^/dz'^ + V{z) with eigenvalues Aq, . . . , A„. Then the differential 
operator 

-^ + V-2(\ogW{i^o,...,^n))", (5) 

where W is the Wronski determinant, has the same eigenvalues as L, except 
Ao, • • • , A„. 
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Let J > be an integer, (&, A) G Z}(R), and let y be the eigenfunction 
corresponding to (6, J, A). Wc apply the Darboux transform to our operator 
Lb^j in (3), taking aU QES eigenfunctions as ipo^ . . . ,i/jn- If J = the Dar- 
boux transform does not change anything. It is easy to see, [12] that the 
transformed operator (5) is in this case. But Lb^j{y) — Ay, and if we 

define y*{z) — y{—z), then L^^_j{y*) — Xy*. However, this y* is not an eigen- 
function of Li, _j, because it does not satisfy the normalization condition in 
(3). Instead it tends to zero in 5*2 and in S'_2. This means that y* is linearly 
independent of the eigenfunction yo of L^^-j, and g = yo/y* has asymptotic 
values oo in S2 and in 5'_2- 

Since g has equal asymptotic values in 5*2 and 5'_2, any other Nevanlinna 
function for the same (6, J, A) has the same property. So c = c at the point 
(6, -J, A). 

The case of negative J is treated similarly, applying the inverse Darboux 
transform. Thus c = c on ^}(R) when J is an integer. This proves the "only 
if" part of Theorem 2. 

As a byproduct we obtain the following: 

(iv) condition c — 00 at the point {b, A) e Z_j(R) holds if and only if J > 
is odd and {b, A) e Z^^^(R) n Z}(R) 

Indeed, the eigenfunction y is elementary if and only if the asymptotic 
values of / in S3,Si,S-i arc zero. This happens if and only if the function g 
defined above has asymptotic values 00 in Sq, S2, S-2- 

To prove the second part of Theorem 2, we need to find which part of the 
spectral locus corresponds to real c. 

4 Line complexes 

Assigning asymptotic values in sectors Sj is not enough to define a Nevan- 
linna function /, one needs additional information about the topology of the 
covering 

/ : C\/~^ (asymptotic values) — >■ C\{asymptotic values}. 

Such information is encoded in the following way [18, 9, 10]. One chooses a 
cell decomposition $ of the Riemann sphere C such that each 2-cell contains 
one asymptotic value, and takes the /-preimage of this cell decomposition. 
This preimage is a cell decomposition ^ of the plane which locally looks 
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like There are many ways to choose $, and here we describe cell de- 
compositions used by Nevanlinna, see also [15, Ch. 7] for a comprehensive 

treatment. 

First we fix two points in C which are distinct from the asymptotic values. 
We call these points x and o. They are the vertices of Suppose that we 
have q asymptotic values. We connect x and o with q edges which do not 
intersect except at the ends. These edges and vertices form the 1-skeleton of 

we require that faces of $ contain one asymptotic value each. This choice 
of $ is fixed in this section. 

The preimage \& = /~^($) is called the line complex. It has the following 
property: the 1-skeleton of is a bipartite connected graph embedded in C 
whose all vertices have degree q, and each component of the complement has 
either two or infinitely many edges on the boundary. Moreover, the number of 
components having infinitely many boundary edges is finite. These properties 
completely characterize all possible line complexes arising from Nevanlinna 
functions. 

This means that for arbitrary line complex, and asymptotic values (Na- 
evanlinna parameters), satisfying the restrictions stated in section 2, there 
exists a normalized Nevanlinna function with this line complex and these 
asymptotic values, and it is unique up to an affine change of the independent 
variable. The Schwarzian derivative of this function is a polynomial and 
we compose with an affine change of the independent variable so that this 
polynomial has leading coefficient —2 and next coefficient vanishing, so that 

^ - 2 (77] = -2(^' + <^<i-2z'-' + ... + ao). (6) 

In this paper, d = A, (02, Oi, ao) = {2b, 2iJ, A) and Nevanlinna parameters 
are asymptotic values a and c as described in section 2. Function / does 
not change if a and c are multiplied by the same number. For a fixed line 
complex, the map that sends the parameter t — c/ato the Nevanlinna triple 
(6, J, A) is called the Nevanlinna map. It is smooth, injective and has non-zero 
derivative. It gives a paramctrization of a part of the spectral locus. 

We label the faces of $ and \I' with the asymptotic values, so that a face 
of ^ has the same label as its image. Two cell decompositions "if' and ^" are 
considered equivalent if one can be mapped onto another by a homeomor- 
phism of the plane preserving orientation and labels of faces and vertices. If 
two Nevanlinna functions / and g have equivalent cell decompositions, then 
f{z)^g{az + /3), a 0. 
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The following properties are evident. The cyclic order of face labels 
around a vertex of ^ is the same as for the image vertex in $. Two faces of 
\E' with the same labels have disjoint closures. Each face is bounded either 
by two edges of by infinitely many edges. 

By erasing multiple edges of the 1-skeleton of ^ and discarding the labels 
of bounded faces, we obtain a new cell decomposition with labeled faces 
and vertices, whose 1-skeleton is a tree. The line complex can be uniquely 
recovered from its associated tree. 

If / is real, so asymptotic values are symmetric with respect to complex 
conjugation, sometimes it is possible to choose a symmetric cell decompo- 
sition $. Then ^ is also symmetric. And conversely: if the labeled cell 
decompositions $ and \& are symmetric with respect to complex conjugation 
then / can be chosen real by pre-composing with an afiine map of C. 

Line complexes are convenient for study of the limits of families of Nevan- 
linna functions when two asymptotic values collide. We have the following 
compactness theorem [22]. Fix a cell decomposition Let be a se- 
quence of Nevanlinna functions with line complexes = /~^($). Suppose 
that Vn = is a vertex of "^n, and that /„ are normalized by conditions 
l/n(0)l — 1- Then one can choose a subsequence from that tends to a 
limit, and this limit is a Nevanlinna function. 

We need only the special case when the cell decompositions corre- 
sponding to /„ are all equivalent. If distinct asymptotic values of /„ tend to 
distinct limits, then / has the same cell decomposition ^. If two asymptotic 
values of fn which are labeling adjacent faces of $, collide in the limit, then 
one has to erase from the 1-skeleton of ^ all edges on the common bound- 
aries of faces with these coUided asymptotic values. The component of the 
remaining graph containing the vertex v is the cell decomposition of the limit 
function. 

Now we return to Nevanlinna functions corresponding to problem (3). 
One technical problem we are facing is that it is not always possible to 
choose a symmetric However this is possible when c is real, and in the 
next sections we consider this case. We begin with the simplest case when 
J = 0. In this case the Nevanlinna function has an additional symmetry. 

The cell decompositions considered in [13] are different from line com- 
plexes, because in the situation considered in that paper, it is impossible to 
define a line complex with the required symmetry properties. 
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5 Subfamily J = 



In this section we begin to prove Theorem 1. 

Let y be an eigenfunction of Lfeo, where 6 G C. Function yi{z) = y{—z) 
satisfies the differential equation in (3) with J = 0, but does not satisfy the 
boundary conditions. So yi is hnearly independent of y, and we consider the 
Nevanhnna function 

/ = y/yi- (7) 

It is not normalized'm the sense of section 2. This function / has the following 
symmetry property: 

/(-;.) = !//(.). (8) 

The asymptotic values are in Si and S_i, oo in 5*2 and S'_2, A in 5*0 and 
1/A in 5*3. This A is the Nevanlinna parameter. It follows from (6) that the 
condition (8) is equivalent to J = 0. 

If h is real, we know from the results of [4] and [20] that all eigenvalues 
of Lq i, are real. Hence / is a real function and A G R. As / is defined up to 
multiplication by a real non-zero number, we can normalize so that A> Q. 

Proposition 1. For 6 G R, we have A G (0, 1). Each eigenfunction has at 
most one zero on the real line. 

Proof. Function / has the property that 

/ : C\/-i({0, A, IM, oo}) ^ C\{0, A IM, oo} 

is a covering. To construct the line complex, we choose the cell decomposition 
$ of the target sphere shown in Fig. 1. 




Fig. 1. Cell decomposition $ of the Riemann sphere. 
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It has two vertices, four edges, and four faces labeled by the asymptotic 
values. We denote by A' the asymptotic value which is in (0, 1), so that A' is 
either A or A~^, and our first goal is to find out which of these possibilities 
holds. 

The line complex ^ = /"^i^) is a labeled cell decomposition of C. It 
has 6 unbounded faces. Moreover, ^ is symmetric with respect to the real 
line and with respect to the imaginary line. The symmetry with respect to 
the real line does not change the labels, while the symmetry with respect 
to the imaginary line interchanges x with o, with oo and A' with 1/A'. 
Unbounded faces of ^ are asymptotic to the sectors Sj. 

For any pair of vertices of ^ connected by several edges, we replace these 
several edges with one edge. The result is a simpler cell decomposition T 
whose 1-skeleton is a tree. The faces of T are labeled with asymptotic values, 
and T has all symmetry properties described above. The label of a face 
asymptotic to Sj is the asymptotic value in Sj. 

It is easy to classify all possible labeled trees satisfying the above con- 
ditions. They all have two vertices of order 4 and the number k of edges 
between these two vertices is odd. These trees depend on one non-negative 
integer parameter m, such that k = 2m + 1, and we denote them by T^. The 
tree Ti and the corresponding line complex are shown in Fig. 2. 
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Fig. 2. The tree Ti (left) and the corresponding hne complex (right). 

Comparing the cyclic order of labels of faces adjacent to a vertex of order 
4 in Fig. 1 and Fig. 2 we conclude that A = A' so A & (0, 1). This proves 
the first part of the proposition. The second part is immediately clear from 
the classification of the trees T^: the cell decomposition \1/ has at most one 
face labeled with which intersects the real line. (A similar argument was 
used in [9].) Notice that 

the number of real zeros of ?/ is if m is even and 1 if m is odd. (9) 

Lemma 1. For b G R, eigenvalues of Lq i, are distinct. 

Proof. The subset of the spectral locus Z parametrized by the Nevan- 
linna functions satisfying (8) coincides with Zq. Since Nevanlinna parameters 
(modulo multiplication by non-zero numbers) serve as a local coordinate sys- 
tem on Z, Zq is a smooth one-dimensional subset of Z. Hence ZqCR,) consists 
of smooth real curves and, possibly, isolated points. Since all eigenvalues of 
h are real for 6 G R, isolated points are not allowed. Hence Zo(R) is 
smooth, and the real eigenvalues of L^fi cannot collide as 6 G R changes. 
This proves the lemma. 
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Thus we can label the eigenvalues in increasing order, Xo{b) < \i{b) < 

To each eigenvalue corresponds one cell decomposition ^, and one tree 
T = Tm- So we have some correspondence m{k), which is defined so that 

Tm{k) corresponds to A^. 

Proposition 2. m{k) = k. So the eigenfunction corresponding to Xk has 
no real zeros when k is even and one real zero when k is odd. 

Proof. We use the asymptotic result in [12, 11] to degenerate (3) as 
h +00 to the harmonic oscillator 

Y"{u) + Au^Y{u) = nY{u), Y{it) ^0, ±oo. 

The boundary condition here comes from the boundary condition in (3). 
Eigenfunctions Yfc corresponding to eigenvalues //fc, labeled in the increasing 
order, have all zeros on the imaginary axis. There is one zero on the real axis 
(namely at the origin) when k is odd and none if k is even. 

To analyse the behavior of zeros of eigenfunctions of (3) as 6 +cx3, 
we consider the function g = f/{f — A). This corresponds to a difi'crcnt 
choice of yi in the basis {y, yi) of solutions of the differential equation in (3). 
Asymptotic values of g are oo in 5*0, in S±i, 1 in S±2, and a :— 1/(1 — A^) 
in 5*3. 

According to the result of [11], to every eigenfunction with eigenvalue 
Ilk, and to every positive b large enough, corresponds a unique eigenfunction 
y of (3) with eigenvalue X{b) — > /ik- Let m — m{k) and let be the tree 
corresponding to y. Then must collapse to the tree corresponding to 
n (see Fig. 3). 
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Fig. 3. Tree T[. 

By collapse of a tree we mean the following: two face labels become 
equal, and all edges on the common boundary of these faces are erased. Such 
collapse can happen only if 1/(1 — A^) — )■ 1 or 1/(1 — A^) — )■ oo. In the second 
case we must multiply by 1 — A^, so the asymptotic value in 5*3 becomes 1 
and in 5*^2 it becomes 1 — — )■ 0. Then as the tree Tm{k) collapses to the 
tree T^, we must have m{k) = k. This proves the proposition. 

Proposition 2 gives a parametrization of the spectral locus for J = by 
Nevanlinna parameter A. For each m, the correspondence 6 — )■ A is a real 
analytic homeomorphism of the real line onto (0, 1). 



6 Region J < 1 

Now we can prove Theorem 1. When b, J and A are real, the eigenvalue 
problem is PT-symmetric, and we can define a real eigenfunction y. We 
choose the second linearly independent solution yi of the differential equation 
in (3) from the conditions that — j- in 5*0, and yi is real. Then f = y/yi is 
real. The asymptotic values of / in 5*0, . . . , 6*5 are cxd, 0, c, a, c, in this order, 
where a G R. Thus for every (6, J, A) on the real spectral locus, Nevanlinna 
parameters a G R and c G C are defined. 
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As long as a 7^ 0, the number of real zeros cannot change, because y and 
/ never have multiple zeros. 

Equality a = can happen on the spectral locus if and only if J is a 
positive integer and the eigenfunction y is elementary [2, 12]. 

Lemma 2. Let j{t),t G [0, 1] be a curve in the real (6, J) plane. Suppose that 
7(0) = (0,0), and that J < 1 on j. Then for {b,J) = 7(1), all eigenvalues 
are real. 

This implies Theorem 1. 

Proof of Lemma 2. It is sufficient to prove the lemma for curves in the 
open half-plane J < 1. The general case follows by continuity. 

For (6, J) = (0, 0) all eigenvalues are real and distinct, so we can order 

them as Ao < Ai < Then each A^ can be analytically continued along 

7 for t e [0,tk) for some tk > 0. We denote these analytic continuations 
by Xk{t)- According to theorems of Shin [19] all but finitely many of the Afc 
are real and by [10, Thm. 1], all but finitely many A^ can be analytically 
continued for t e [0, 1]. Let G be a bounded simply connected neighborhood 
of 7 in the half-plane {{b, J) E : J < 1}. Then all but finitely many 
branches of X{b, J) are holomorphic, distinct and real in G. The remaining 
branches satisfy a minimal algebraic equation of the form 

A"^ + am-i{b, J)X^-^ + . . . + a„(6, J) = 0, (10) 

with aj analytic in G. The zeros of the discriminant of this equation in G 
form a closed subset K d G. 

We claim that all solutions of equation (10) are distinct for all (6, J) e 7. 
Indeed, if any two eigenvalues collide as t increases, then some adjacent eigen- 
values \j{t) and Aj+i(t) must collide at some point to- The corresponding 
eigcnfunctions will tend to the same limit as t t^. But this is impossible 
because one of them has no real zeros and another has one. 

Thus K does not intersect 7, and there is an analytic continuation of 
eigcnfunctions to 7(1). 

7 Classification of line complexes 

In this section we consider the general case of real asymptotic values a and 
c. In this case the normalized Nevanlinna function / has asymptotic values 
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(oo, 0, c, a, c, 0) in 5*0, ... , S^. We choose the cell decomposition $ similar to 
that in Fig. 1, see Figs. 5,7,9. 
There are three generic cases: 

Case =Sf . c < < a. Case < a < c, Case tf. < c < a. 

Letters J^, ^ and S stand for "left", "right" and "even", the meaning of 
this notation will be clear later (Figs. 11-15). There are also non-generic 
cases a = and c = oo. Assuming that a 7^ 0, we classify all possible line 
complexes. 

In all cases =Sf, =^ and S ^ we first classify all possible bipartite trees 
symmetric with respect to the real line, with 6 faces labeled by 00, 0, c, a, c, 
in this cyclic order, the face labeled 00 bisected by the positive ray, and 
satisfying the condition that faces with the same label have disjoint closures. 
There are three types of such trees shown in Fig. 4. They depend on two 
integer parameters. A; > and /, where k is the number of edges between two 
ramified vertices as shown in Fig. 4. Parameter I takes all integer values, but 
|Z| is the number of edges between ramified vertices as indicated in Fig. 4. 
We say that a tree has type if Z = 0, type 1 if / > and type 2 if Z < 0. 

For the trees that occured in section 5 we have T„ = X2„+i,o- 



/>0 







/=0 



c 




00 



Fig. 4. Three types of trees. 
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Now we consider all cases separately, and argue by the following scheme. 
First, for a given case of ordering a,c, on the real line, and for each tree, 
we decide whether this tree can come from a line complex, and if it can, we 
recover the line complex. 

To do this, we begin with a ramified vertex v of the tree. Comparing the 
cyclic order of the labels around this vertex v with the cyclic order around 
the vertices of the cell decomposition $ of the sphere we determine whether 
this vertex of the tree is a x or o, and add the missing edges, step by step. 
When we come to another ramified vertex, either the cyclic order is correct 
or not. If it is not correct, the tree does not correspond to a line complex in 
the considered case. Otherwise, we recover the line complex uniquely. 

We write the Nevanlinna map defined in section 4 as (6, J, A) = 
where ^ is the line complex, and t = c/a. 

Then we consider possible limits as a — or c — )■ oo on the spectral locus. 
This gives degenerations to the intersections of the QES spectral locus with 
the non-QES spectral locus, and description of these intersections in [13] 
permits to recover the value of J from the tree in cases 1 and 2. Then we 
consider the degeneration as c — )■ oo or c — )■ a whenever possible, as it was 
done in [11, 12, 13]. We will conclude that real values of c correspond to 
integer J and we obtain the parametrization of the whole spectral locus for 
integer J. 

Case =Sf . c < < a, see Figs. 5-6. 




Fig. 5. Cell decomposition $ in case =Sf. 
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Fig. 6. S^m.i complexes. 

Trees of type (with / = 0) are impossible in this case. 

Trees X^, ; of type 1 (with / > 0) are possible in this case if and only if 
k = 2m + 1, m = 0, 1, . . . and / > 1 is odd. Such line complexes will be 
called ^m,i- When a \ 0, t = c/a — )■ — oo, the corresponding Nevanlinna 
function has a limit on the QES locus. This limit function has / — 1 zeros, 
none of them real. It follows that J = I for the limit function. By the first 
part of Theorem 2 which we proved in section 5, c is real for integer J on the 
non-QES part of the spectral locus, so we conclude that the whole image of 
the Nevanlinna map t t— )■ F{^rn,iit) belongs to ^j(R) with J = I. 

In the limit when c 0, that is t = c/a 0, we obtain a Nevanlinna 
function for the harmonic oscillator with m zeros. 

Trees of type 2 (with / < 0) are possible in this case if and only ii k = 
2m + 1 and / < — 1 is odd. Such line complexes will be called ^m,i- The 
only possible limit on the spectral locus is c — )■ — oo. This corresponds to an 
elementary second solution of the differential equation in (3) (solution which 
is linearly independent of the eigenfunction). These points are marked in 
Fig. 15. Thus the chart ^m,i, ^ < — 1 corresponds to the chart S^m-i via 
Darboux transform. We have J = / < in this case. 



Case ^. < a < c. Figs. 1-i 
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Fig. 7. Cell decomposition $ in case ^. 
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Fig. 8. <^m,z complexes. 



Trees of type (with Z = 0) are impossible in this case. 

Trees ^ of the type 1 (with I > 0) are possible with k = 2m + 1 
and / positive odd. We call the complex ^m,i- Degeneration a — )■ 0, t = 
c/a ^ +00 is possible, and the limit belongs to the QES spectral locus. The 
limit Nevanlinna function has / — 1 zeros, so J = I for this limit function. 
Again, using the first part of Theorem 2, we conclude that J = / on the 
whole image of the Nevanlinna map F{Mm,i, ■) Degeneration c — ?■ +oo (that 
is t = c/a — 7- 0+) gives a Nevanlinna function for the harmonic oscillator 
with J — 1 + m zeros. As =Sfm,z and ^rn,i have common limit of the QES 
spectral locus, their Nevanlinna images form a single curve. The results in 
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[13] about QES spectral locus together with counting of zeros of degeneration 
to harmonic oscillator show that -f^m^i lies on the left and ^m,i lies on the 
right from the intersection point with the QES spectral locus. See Figs. 11, 
13, where Nevanlinna images of =Sfm,« and ^ra,i are shown with the solid lines, 
and the QES locus with the dotted line. 

Trees ^ of type 2 (with / < 0) give line complexes when k = 2m + 1 
and / negative odd. We call these complexes ^rn,i- Degeneration c — )■ oo is 
possible on the spectral locus. These charts correspond to the charts ^m,-i 
by Darboux transform. 

Thus cases =Sf and ^ and trees of types 1 and 2 cover all cases when J is 
odd. We conclude from this: a) J is constant on the Nevanlinna images of 
^m,i and ^rn,i, namely J = I. b) Even values of J must be covered by the 
remaining trees from our classification. 

Case S . < c < a, see Figs. 9-10. 
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Fig. 9. Cell decomposition $ in case S. 
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Fig. 10. I complexes. 



^x 



Trees X^ q of type have parameter k = 2m+l, m > 0, and we denote the 
corresponding hne complex by ^m,o- This line complex represents Nevanlinna 
functions from section 5 which correspond to J = 0. 

Trees Xk^i of type 2 have k = 2m + 1 and I negative even. We call 
the corresponding line complex (S^^i- No degeneration on the spectral locus 
is possible. We know that for even J the non-QES spectral locus consists 
of graphs of functions. Degeneration as c — )■ and c — )■ a gives a Nevan- 
linna function for the harmonic oscillator with m zeros. Thus these trees 
parametrize the whole spectrum for negative even J. 

Trees ; of type 1 have k = 2m + 1 and / positive even. We call 
the corresponding complex S'^^i. It corresponds to <%n,-/ by the Darboux 
transform. Degeneration as c — ?■ and c — >■ a gives a Nevanlinna function 
for the harmonic oscillator with m and m + / zeros, respectively. 

These arguments show that real c and a ^ correspond to integer J and 
that we obtain a parametrization of the whole non-QES spectral locus in this 
way. This completes the proof of Theorem 2. 

The parametrization of the spectral locus for integer J is represented in 
Figs. 11-15. The symbols of line complexes are written below the corre- 
sponding curves. The QES spectral locus is shown with dotted lines. It was 
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parametrized with different cell decompositions (not with line complexes!) 
in [13]. Symbols X^^i in the figures refer to the charts on the QES locus 
described in [13]. 





Fig. 14. Zo(R). 
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^-1,0 



Fig. 15. Z_i{R) superimposed with (R) (thin dotted hne). 
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